Photovoltaic Current Response of Mesoscopic Conductors to Quantized Cavity Modes 
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We extend the analysis of the effects of electromagnetic (EM) fields on mesoscopic conductors to include the 
effects of field quantization, motivated by recent experiments on circuit QED. We show that in general there is 
a photovoltaic (PV) current induced by quantized cavity modes at zero bias across the conductor. This current 
depends on the average photon occupation number and vanishes identically when it is equal to the average 
number of thermal electron-hole pairs. We analyze in detail the case of a chaotic quantum dot at temperature T s 
in contact with a thermal EM field at temperature Tf, calculating the RMS size of the PV current as a function 
of the temperature difference, finding an effect ~ pA. 
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Many quantum electronic devices for applications in 
metrology and quantum information technology involve the 
interaction of electrons with high frequency electromagnetic 
(EM) fields, often the quantum devices act as detectors of this 
radiation d. I n phase-coherent (mesoscopic) devices there 
are quantum interference effects in electron transport such as 
the weak localization correction to the conductance and uni- 
versal conductance fluctuations which can in principle be used 
to detect radiation since it suppresses these effects y|,|3|,|4|]. In 
practice the suppression of coherent transport by EM fields 
is difficult to separate from the suppression by intrinsic inter- 
actions due to the electron-electron and electron-phonon cou- 
plings. 

A more reliable mean of using mesoscopic conductors to 
detect EM radiation is to look at the DC current induced 
by such a field at zero voltage and temperature bias across 
the device, known as the mesoscopic photovoltaic (PV) ef- 
fect 0. 01- This effect arises in mesoscopic conduc- 
tors because the phase-coherent transmission through the de- 
vice almost always violates parity symmetry and the non- 
equilibrium distribution created by the EM field sets up a 
steady-state current dictated by this parity violation. When 
the parity violation is due to random interference, the sign of 
this current will fluctuate from sample to sample and its root- 
mean-square (RMS) size in this case depends on the power in 

the EM field lUElliUHEIl- Hence after this PV current 
is calibrated it can be used for detection of the power in the 

incident EM field. 

The previous theoretical description S SHE EE! 
of the PV current has employed a classical treatment of the 
EM fields, since this description was sufficient for the sys- 
tems studied experimentally . In this case the 
RMS PV current is a monotonically increasing function of the 
EM field power. Recently a new generation of electronic cir- 
cuits was developed 1 16], where a quantum electronic device 
is coupled to an EM field of a high quality electromagnetic 
resonator. If the resonator contains a small number of photons 
of the EM field and the lifetime of the photons is long, the 
interaction of the EM field with electrons requires a full quan- 
tum treatment, based on the laws of quantum electrodynam- 
ics, leading to a new sub-field of quantum electronics known 
as circuit QED. 




FIG. 1 : a) Schematic dependence of the PV current of a mesoscopic 
conductor on the average photon number N f in a single mode res- 
onator. Narrow lines represent magnitude of the current for various 
conductor realizations. The current vanishes for all realizations for 
Nt = NsiLjj/T). Bold line is the RMS value of the current over 
the ensemble, calculated for the case of a chaotic quantum dot. b) 
A possible experimental setup for the observation of the quantum 
photovoltaic effect. Two gates of a quantum dot are terminals of a 
high quality microwave line of length L, that form an EM resonator. 
To reduce dissipation, the substrate of the resonator contains 2DEG 
(brighter area) only in the vicinity of a quantum dot. 



In this Letter we investigate the properties of the meso- 
scopic PV current that arises due to the electron interaction 
with quantized EM fields. We find that the net current through 
the device can be represented as the sum of two opposite 
contributions. One contribution is determined by the aver- 
age number N, of photons with energy ai,- in modes i of the 
resonator. The second contribution is determined by the num- 
ber of thermal electron-hole pairs with energy co, of resonator 
modes i, given by the Bose distribution function NB((L>i/T e ) at 
temperature T e of electrons in the leads; Nb(x) = l/(e' 1 ' - 1). 
We demonstrate that if both contributions are taken into ac- 
count, the magnitude of the PV current, unlike the classical 
case, is not monotonic function of the strength of the EM 
field. Instead special conditions can be met when these two 
contributions cancel each other and the PV current vanishes 
for all mesoscopic realizations of the device. For the case 
of a thermal photon field the zero current state occurs when 
the temperature of the EM field and that of the electrons in 
the leads are the same and follows from the principle of de- 
tailed balance. For an externally driven single-mode cavity 
the zero current state also occurs whenever the average num- 
ber of photons Nj is equal to the occupation number of bosons 
A^B(w,/r e ) at the electron temperature T e , independent of the 
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other properties of the full photon distribution in the cavity. 
In the classical description of EM fields only the first contri- 
bution is found and the second contribution due to electron- 
hole pairs is missed. As noted, in this case the PV current 
never vanishes simultaneously for all realizations and its RMS 
value is a monotonic function of the power of the EM radia- 
tion 1 51 la, UjJ, 1 1 '711 . The non-monotonic behavior of the RMS 
value of the PV current is an indication of the quantum behav- 
ior of the EM fields. 

A schematic depiction of the dependence of the PV cur- 
rent on the number iV; of photons in a single mode resonator 
is shown in Fig.^^ f° r the case of any mesoscopic conduc- 
tor with fluctuating transmission matrix. For a particular re- 
alization of the conductor, the current is a linear function of 
Ni and changes its direction at N, = NB(cOi/T); at large val- 
ues of N( the dependence on jV, will depart from linearity due 
to the suppression of the coherence time by electron-photon 
scattering. We note that the average value of the current with 
respect to realizations of the mesoscopic conductor is zero, 
since the parity violation of transmission is zero on average. 
In this case we characterize the magnitude of this current by 
the RMS value of this current averaged over the realizations, 
shown in Fig.^a by a bold line. 

The generation of the photovoltaic current, studied here, is a 
common phenomenon for out-of-equilibrium mesoscopic sys- 
tems, which function as "quantum ratchets". Other similar 
phenomena are the Coulomb drag current fljll and the cur- 
rent in mesoscopic metal rings coupled to out-of-equilibrium 
electron 11911 orphonon 1 20] reservoirs. 

Model — The signatures of the quantum behavior of the 
EM field can be observed in the PV current measurements for 
various mesoscopic systems, such as quantum point contacts 
0], metal ringsjjl2il. metal wires or grains, and semiconductor 
quantum dots |6]. The main requirement is that the device has 
long coherence and inelastic relaxation times, so that electron 
interference in propagation through the device leads to strong 
energy dependence and intrinsic parity violation in transmis- 
sion to the left and right lead. Any such mesoscopic device 
will show a PV current with the properties depicted in Fig. 
□a. 

The specific case we will now treat in detail is the PV cur- 
rent through a semiconductor quantum dot with a few open 
channels placed at the terminal of an electromagnetic res- 
onator, see Fig. ^ b. The quantum dot is similar to those 
used in charge pumping experiments H5T1 . The electro- 
magnetic resonator consists of a microwave line of length L, 
characterized by high quality factor and resonant frequencies 
a>i = inc*/L 11611 . where c, is the phase velocity of EM wave 

in the resonator and i = 1,2, The Hamiltonian of this 

system is 



<H = <H d + <Hi + % + Hd, 
where ^Ha is the Hamiltonian of the electrons in the dot 



n,m= 1 



(1) 



(2) 



Here t]/ n and a, are the annihilation operators of electrons in 
the dot in state n (n — 1 , . . . , M) and photons in mode i of the 
EM field, M x M Hermitian matrices H and V, represent the 
stationary part of the electron Hamiltonian and the electron 
coupling to mode i of the EM field, respectively. 'Hf describes 
the evolution of the electro-magnetic field and can be written 
in terms of photon annihilation and creation operators a, and 



(3) 



where w, is the energy of photon excitations. 

The Hamiltonian for electrons in the leads near the Fermi 
surface is 



(4) 



where i// a (k) is the annihilation operator of electrons in chan- 
nel a of one of the leads. The continuous variable k denotes 
electron momenta in the leads, up = (2ttv) _1 is the Fermi ve- 
locity, and v is the density of states per channel per spin at the 
Fermi surface. In this Letter we consider the case when the 
voltage bias across the dot is zero. The coupling of electron 
states in the dot to states in the leads can be written as 



Hi 



= Y J (w na ^t(k)^ n + n.c). 

a, n.k 



(5) 



Here a labels channels in the leads, with 1 < or < N\ for the 
iVi channels in the left lead and with N\ + 1 < a < N c h for the 
Nr channels in the right lead, A^h = N\ + N T . The coupling 
between electron states in the leads and in the dot is described 
by A^ c h x M matrix W. 

Photovoltaic current (PV) — We calculate the PV current 
that flows through a quantum dot at zero temperature and 
voltage biases. The interaction of electrons with the EM 
field results in the deviation of the electron distribution func- 
tion «d(e) in the dot from the Fermi distribution function 
«f(£) = [1 + exp(e/r e )]~ 1 of electrons in the leads at temper- 
ature T e and in a finite electric current through the quantum 
dot. The direction and the magnitude of such current depend 
on the mesoscopic violation of the left-right symmetry of the 
dot, on the electron spectrum in the dot, and on the coupling 
strength of electrons to the EM field. The derivation of the 
expression for the current follows along the lines for the cal- 
culation of the current through open quantum dots coupled to 
classical external fields |6]. In the case of quantum fields, the 
field acquires the off-diagonal matrix elements in the Keldysh 
space, which can be easily taken into account within a bilinear 
response. As a result, we have 



I = e^^ J Ji(s, ±iodRi{s, ±iOi)ds 



(6) 



The kernel J,(e, oj) contains all the information about electron 
motion in the dot 

JiQs, u>) = 4-7tv ti[WAW t G r (e)y/Im{G r (e - w)}V,G a (e)] (7) 
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and corresponds to the triangle vertex diagram for the 
Coulomb drag 1 18], written for the open dot geometry. Func- 
tion G T (s) is defined for a given realization of H by 



G r (e) 



1 



s-H- invwW* 



G.(e) = [G r (£)f 



(8) 



Here A = (N r /N^ h )Ai-(Ni/N ch )A r , where [A,]^ = 5 afj for 1 < 
a,B < N\, and [Ai] a « = otherwise; A r = 1 - K\. Equation l|6} 
takes into account spin degeneracy. 

The function Ri(e, a>i) is a combination of the Fermi n-p and 
Bose Nb functions, and the photon occupation number A',- = 
(a* a,) for the i mode of the EM field: 

Ri(s, u) = 2 [Nt - N B (a>/T e )] [n F (e - u) - n F (e)] . (9) 

According to Eq. the contribution to the current / from 
mode i of the EM field is linear in the average number A, of 
photons with energy w, in mode i. In particular, if only one 
mode of electromagnetic field is coupled to electrons in the 
dot, the current / can be used to determine the average number 
of photons in this mode. 

The structure of Eq. can be understood from the follow- 
ing schematic argument. The distribution function in the dot 
«d(e) is the solution of the kinetic equation 



«dO) - «fO) 



2[rf( e )-rn £ )]. (io) 



Here the left hand side describes the relaxation of the distri- 
bution function «d(s) due to electron escape to the leads with 
characteristic escape time T esc and the right hand side repre- 
sents the imbalance between the rates of absorption and emis- 
sion of photons. These rates are determined by «d(e) and by 
the average number of photons A, in mode i of the electro- 
magnetic field: 



TfocNmieXl-mie + (*)), 
IT « [#+l]Bd(e + «■*)(! -nd(e)). 



(11) 



Equations d 1 0t and dl 11 1 determine the distribution function 
in the dot «d(e). To the lowest order in electron coupling to 
the EM field, we can substitute «d(e) = rip(s) in Eq. dl 11 1 and 
obtain [r*(e) - r™(s)] oc 7?,(e, w,). A similar structure of the 
expression for the electric current was also obtained in 12011 
for electron system coupled to an out-of-equilibrium phonon 
reservoir. 

When the thermal state of the EM field at temperature Tf 
is equal to electron temperature T e , A, = AB(w//T e X the so- 
lution of Eq. fllUi is «d(e) = n F (s). The latter equality is the 
consequence of the detailed balance principle. Since the dis- 
tribution functions in the dot and in the leads are equal, the 
PV current vanishes, as expected for a system in full thermo- 
dynamic equilibrium. 

Note however that the PV current does not vanish when 
there are no photons in the cavity, A, = 0. In this case the 
current is driven by electron relaxation through spontaneous 




FIG. 2: False color plot of the kernel "K determining the RMS PV 
current as a function of its arguments, the cavity mode frequency, 
a>[, and the electron temperature, T c , in units of the electron escape 
rate y esc from the dot. The maximum RMS current will occur at the 
maximum of "K, which occurs when <±>i,T c ~ y esc . 



emission of photons to the unoccupied modes of the EM field, 
and the magnitude of this current is determined by the matrix 
elements V, of zero-point fluctuations of the EM field. 

Below we consider the case when the electromagnetic field 
is in a thermal state at temperature Tf and A, = Nn{o)ilTf). 
Function Ri(e,co) contains Nn(coilTf) - AB(w ; /r e ) and van- 
ishes identically for Tf = T e . For small deviations of Tf 
from T e , the current / is linear in the temperatures difference: 
/ = B(Tf - r e ), where the EM field thermopower coefficient B 



B 



r(s, oS) — 



(s, Lu>i)r(s, (of) - Ji{s, -(L>i)r(s, -o»,)] da 
w/[4r e 2 sinh(w/2r e )] 



cosh(e/2r e ) cosh[(e - (L>)/2T e 



(12) 



Equation dl 2b determines the value of the EM field ther- 
mopower coefficient B for a particular realization of the 
Hamiltonian H of the quantum dot. Theoretical and exper- 
imental work has shown that lateral quantum dots 
are well-described by a random-matrix (RM) model of their 
transmission properties due to the chaotic motion of electrons. 
Therefore the theory of such systems has focused on calculat- 
ing averages of relevant statistical quantities over an appropri- 
ate RM ensemble. Below we calculate the statistical proper- 
ties of B with respect to a RM ensemble of H. 

Mesoscopic fluctuations of the current — We calculate 
mesoscopic fluctuations of thermopower coefficient B with re- 
spect to realizations of the M x M matrix H from a Gaussian 
ensemble of Hermitian matrices with M — > oo, characterized 
by the mean level spacing 5,: (H mn H* nm ) = (M5\ln 2 )5 nnl 5 mm , 
for a unitary ensemble [The result is the same for unitary and 
orthogonal ensembles]. The quantity of interest is the dimen- 
sionless quantity representing the ensemble average product 
of interaction matrices V,: 



T ij = 



It: 1 <tr{fcfy}> e 



(13) 



A ch (Mdi) 2 

The coupling constants between electron states in the dot and 
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in the leads are W„ a = 6„ a ^M5\ln 2 v and result in the fol- 
lowing value of the electron escape rate from the dot: y esc = 
1 /i"esc = A^ c h(5i jln. Note that y esc plays the role for this ballis- 
tic system of the Thouless energy of diffusive systems, it sets 
the scale of variation of the transmission with electron energy. 

The ensemble average value of B is zero. The variance of 
the PV current can be calculated by diagrammatic RM tech- 
nique 1 23] and is given by 



2 aw v r„r 

var B - e — ) 

ch ij 



J.I , 



Tesc 



(14) 



where the kernel 7C /; is given by 



Kco,,coj =4y; 



3 + 



2[r?,/(r !i r^.)]^ sc 



yl sc + (S - S' + LOj - LOj) 2 

r{e, a>i)r(s, uj) 



(15) 



y 2 sc + (e - e') 2 



dede . 



To analyze the properties of the variance of the thermo- 
electric coefficient B, we consider the case when only one 
mode i = 1 of the electromagnetic field is coupled to elec- 
trons in the quantum dot. We describe the properties of 
"7Cn = 'K{co\ljesc, r e /y esc ), the contour plot of this function 
is shown in Fig. 13 

In the low frequency limit of co\ <K y esc , the power law 
7C(a»i/y e sc, T e /y esc ) ~ co\ is similar to the dependence of the 
variance of the PV current induced by a single-parameter clas- 
sical perturbation |6]. At low temperature T e <k a>\, the 
number of photons and electron-hole pairs is exponentially 
suppressed and 7C oc exp(-«i/T e ). At high temperature 
T e » y esc , the contribution to the thermoelectric coefficient 
comes from electron states within thermal energies and be- 
comes self-averaged. As a result of such self-averaging, the 
variance of B decreases as 7C ~ 1 /T e as T e increases. To sum- 
marize, < K(u>\ /y e sc, r e /y eS c) has a maximum at T oc u\ at fixed 
(Oi, see Fig. |2 The global maximum of < K{(i)\ly e so, T e /y eS c) 
is 7C max ~ 12.5 at oj\ * 1.2y esc and T e as 0.6y esc - Thus, the 
largest effect will be observed when T e , co\ and y esc are all of 
the same order of magnitude. 

Conclusions — We discuss experimentally achievable val- 
ues of the system parameters (restoring ft and k% in the equa- 
tions below). In experiments 11611 to\ /2n ~ 10GHz (hcai/ks ~ 
0.5K) and 7f w 30mK. The escape rate y esc for N c h ~ 1 is 
comparable with 6i/2nti * 2.5GHz (5i/k B * 0.12K) Q. 

To estimate the RMS value of the thermoelectric coefficient 
firms = VvarB, we write B lms ~ eT\i(k^ln), see Eq. dl4t . 
where T\\ is defined by Eq. Jl 3i and can be expressed in 
terms of the magnitude of zero-point electric field E\ of the 
lowest frequency mode i = 1 of the EM resonator as Fn * 
e 2 E 2 a 2 T l T e s C /ti 1 1 23], where a is the diameter of the dot, see 
Fig. n an d T t = a/vp is the traversal time. The zero-point 
electric field E\ can be estimated from the following equa- 
tion (E 2 /4n)d 2 L ~ h(L>\/2, where L = nc*/u>i is the resonator 



length, and gives the estimate B rms ~ ea t (ks/h)(a> 2 T t T esc ), 
where a* = e 2 /hc, and a ~ d. Compared to the usual 
thermopower due to the temperature difference between the 
leads 1 24], the EM thermopower at T e =s Hy esc /kB is sup- 
pressed by factor a t T t y esc <C 1. At cj\l2n = 10GHz and r t = 
4 • l(T 12 s, we have B rms ~ 35pA/K, and for \T f - T e \ * O.IK 
the current ~pA is in the observable range. 
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